ELSEVIER 


Available  online  at  www.sciencedirect.com 


SCIENCE 


DIRECT® 


Journal  of  Power  Sources  158  (2006)  663-672 


www.elsevier.com /locate /jpowsour 


Predicting  current  flow  in  spiral  wound  cell  geometries 

Jan  N.  Reimers* 

E-One  Moli  Energy  ( Canada )  Ltd.,  R&D,  20  000  Stewart  Crescent,  Maple  Ridge,  BC,  Canada  V2X  9E7 

Received  14  July  2005;  received  in  revised  form  2  August  2005;  accepted  3  August  2005 

Available  online  25  October  2005 


Abstract 

A  general  method  is  described  for  calculating  potential,  current  and  heat-power  distributions  in  electrode  foils.  The  resulting  theory  is  suitable 
for  cell  designs  with  an  arbitrary  number  of  current  collecting  tabs.  Example  calculation  results  for  a  three  tab  cell  design  with  an  ohmic  stack 
impedance  are  shown.  Internal  mixing  currents,  a  non-uniform  stack  current  distribution,  and  sharply  peaked  heat  generation  near  the  tabs,  are 
predicted. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Most  commercial  Li-ion  cells  consist  of  layers  of  foil  current 
collectors,  active  electrodes,  and  separator  coiled  into  a  spiral 
geometry,  which  is  often  referred  to  as  a  jelly  roll.  The  jelly 
roll  cell  geometry  is  not  unique  to  the  Li-ion  chemistry.  Emerg¬ 
ing  from  a  jelly  roll  there  must  be  at  least  two  tabs  (usually 
metal  strips)  which  are,  unfortunately,  often  also  referred  to  as 
current  collectors.  During  discharge  of  such  a  cell,  electrical  cur¬ 
rent  will  flow  into  the  anode  tab(s),  along  the  anode  foil,  across 
the  electrode  stack,  along  the  cathode  foil,  and  out  the  cathode 
tab(s). 

Over  the  past  20  years  great  progress  has  been  achieved  in 
understanding  the  electrical  and  ionic  current  flows  across  the 
electrode  stack  of  Li-ion,  NiMH  or  lead  acid  cells.  Given  a  suffi¬ 
cient  quality  of  transport  and  thermodynamic  properties  for  the 
electrodes  and  electrolyte,  it  is  possible  to  predict  the  dynamic 
response  of  the  stack  for  a  given  stack  design.  Having  this  capa¬ 
bility  should  be  an  integral  part  of  any  cell  design  engineer’s  tool 
kit.  The  theory  and  practical  implementation  of  these  techniques 
has  been  described  elsewhere  (for  a  comprehensive  review  see 
ref.  [1]  and  references  therein). 

Much  of  this  work  treats  the  electrode-separator  stack  using  a 
one  dimensional  model,  the  assumption  being  that  the  electronic 
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conductivity  of  the  current  collector  foils  is  very  high  in  relation 
to  the  effective  conductivity  of  the  stack.  Hence  any  potential 
drop  in  the  foils  and  tabs  can  be  ignored.  This  is  indeed  reason¬ 
able  for  many  commercial  cell  designs.  However,  with  recent 
interest  in  power  optimized  cells  aimed  at,  for  example,  HEV 
applications,  the  geometry  may  be  modified  to  the  point  where 
the  impedance  of  foils  becomes  significant. 

Some  previous  workers  in  this  field  have  addressed  this  issue 
in  their  models  [2-4]  for  simple  tab  arrangements  (one  anode 
and  one  cathode  tab).  The  purpose  of  this  work  is  to  outline  a 
practical  method  for  calculation  of  foil  current  and  potential  dis¬ 
tributions  for  the  general  case,  with  any  number  of  tabs,  located 
at  arbitrary  positions  along  the  length  of  the  electrode  foils.  With 
commercial  cells  already  on  the  market  using  numerous  tabs  per 
electrode,  this  is  not  just  a  purely  academic  issue. 

2.  Cell  geometry 

Unfortunately  the  term  current  collector  is  often  overloaded 
to  refer  to  both  the  tabs  and  the  electrode  foils.  Since  the  dis¬ 
tinction  is  important  for  understanding  this  work  usage  of  this 
overloaded  term  will  be  avoided.  The  cell  geometry  considered 
in  this  work  is  shown  schematically  in  Fig.  1  for  a  non  trivial 
arrangement  of  the  tabs  (two  anode  and  one  cathode  tab).  With 
the  z  direction  going  into  the  page,  the  stack  layers  (foils,  elec¬ 
trodes  and  separator)  are  in  the  xz  plane.  This  diagram  represents 
an  jelly  roll  as  viewed  end  on  after  it  has  been  unwound.  The  y 
direction  has  been  greatly  magnified  for  clarity. 
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Nomenclature 

a ,  b  adjustable  parameters  in  the  general  solution  to  the 

foil  potential  field  equation,  determined  by  exter¬ 
nal  boundary  conditions 

c  adjustable  parameter  in  the  general  solution  to  the 

foil  potential  field  equation,  determined  by  exter¬ 
nal  boundary  conditions  (V  cm-1) 
d  adjustable  parameter  in  the  general  solution  to  the 

foil  potential  field  equation,  determined  by  exter¬ 
nal  boundary  conditions  (V) 
i  stack  current  density  (mA  cm-2) 

I  externally  applied  current  (A) 

I  a,  Ic  anode  and  cathode  foil  currents  (A) 

L  length  of  the  electrodes  in  the  v  direction  (cm) 

n  total  number  of  tabs  sticking  out  of  the  jelly  roll 

n\,  nlc  number  of  internal  tabs  on  the  anode  and  cathode 
electrodes 

neA,  nQc  number  of  external  tabs  on  the  anode  and  cathode 
electrodes 

ha,  nc  number  of  tabs  (internal  and  external)  on  the  an¬ 
ode  and  cathode  electrodes 

n1,  nQ  number  of  internal  and  external  tabs  on  all  elec¬ 
trodes 

0()  higher  order  terms 

q  state  of  charge  in  the  active  stack  (A  h  cm-2) 

qA ,  qc  heat  generation  rate  in  the  anode  and  cathode  foils 
(Worn'2) 

^ stack  heat  generation  rate  in  the  stack  (W  cm-2) 
rA,  rc  relative  foil  impedances,  ta  +  rc  =  1 
Ra  ,  Rc  length  specific  anode  and  cathode  foil  resistivities 
(mQ  cm-1) 

externally  observed,  finite  difference  cell 
impedance  (m£2) 

^ceii  externally  observed,  differential  cell  impedance 

(m£2) 

7? stack  effective  stack  impedance  (£2  cm2) 

t  time  (s) 

7a,  7c  anode  and  cathode  foil  thicknesses  (pan) 

Vceii  externally  observed  cell  voltage  (V) 

V0(q)  open  circuit  voltage  curve  for  the  active  stack  (V) 
Wa,  Wq  anode  and  cathode  foil  widths  (cm) 


( Pa ,  0c  anode  and  cathode  foil  potentials  (V) 

0P  particular  solution  to  the  foil  potential  field  equa¬ 
tion  (V) 

0left ,  0nght  solutions  to  the  homogenous  foil  potential  field 
equation  (V) 

Pa,  Pc  specific  resistivity  of  the  anode  and  cathode  elec¬ 
trode  foils  (jx£2  cm) 

a  electrical  conductivity  in  the  solid  phase  of  an 

electrode  (S  cm-1) 

Subscripts 
A  anode 

C  cathode 

j  segment  or  tab  index 


The  relevance  of  the  jelly  roll  geometry  to  this  work  lies 
only  in  the  length  of  the  electrodes  in  the  v  direction.  Matters  of 
curvature  are  not  important  for  the  calculation  of  dc  foil  current 
distributions  (the  dc  resistance  of  a  piece  of  wire  does  not  change 
when  coiled). 

Current  runs  along  the  foils  in  the  x,  or  in-plane ,  direction  and 
through  the  stack  in  the  perpendicular  y,  or  through-plane,  di¬ 
rection.  Following  the  convention  in  the  cell  modeling  literature, 
a  positive  stack  current  flows  in  the  same  direction  as  positive 
ions,  from  anode  to  cathode,  as  indicated  by  the  vertical  arrows, 
and  corresponds  the  discharge  of  the  cell.  Along  the  foils,  pos¬ 
itive  current  flows  in  the  positive  v  direction,  as  indicated  by 
the  horizontal  arrows.  With  these  conventions,  electron  flow  in 
the  foils  is  of  course  opposite  to  the  current  flow.  In  a  real  cell 
the  extent  of  the  tabs  in  the  v  direction  (tab  width)  is  usually 
<1%  of  the  total  electrode  length.  Indeed  it  is  the  width  of  the 
contact  points  (weld  pattern)  that  is  relevant,  and  this  is  often 
significantly  narrower  than  the  tab  itself.  In  order  to  simplify  the 
analysis  the  contact  width  of  the  tabs  is  considered  to  be  zero. 
The  z  direction,  going  into  the  page,  corresponds  to  the  width  of 
the  electrodes  or  the  length/height  of  the  jelly  roll.  Currents  and 
potentials  are  assumed  to  be  constant  in  the  z  direction.  Again 
this  is  a  very  good  approximation  as  long  as  the  tab  resistance 
in  negligible  in  the  z  direction. 


v  distance  along  the  length  of  the  electrodes  (cm) 

y  distance  across  the  stack  (pan) 

z  distance  across  the  width  of  the  electrodes,  or  jelly 

roll  height  (cm) 

Greek  letters 

a  inverse  length  scale  (cm-1) 

6s°iid ,  6pore  solid  and  pore  volume  fractions  in  an  electrode 
or  separator 

k  electrolyte  conductivity  (S  cm-1) 

0  potential  difference  across  the  stack  under  load 

(V) 


Cathode 


I  >0  Tab  /,.<  0 

c  A  C 


X 


Fig.  1.  Schematic  cell  geometry  corresponding  to  an  unwound  jelly  roll.  Black 
rectangles  attached  to  the  foils  indicate  the  position  of  the  tabs. 
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Stack  #1 


Stack  #2 


Anode 
Separator  #1 

Cathode 

Separator  #2 
Anode 


Fig.  2.  Current  flow  geometry  for  a  spiral  wound  jelly  roll  with  double  sided 
electrodes. 


Fig.  1  is  a  simplified  version  of  the  real  jelly-roll  geometry. 
In  commercial  cells  both  electrodes  are  coated  on  both  sides, 
which  means  each  foil  has  two  stack  directions  to  send  current 
into.  The  cell  actually  has  two  separate  stacks  in  parallel.  Fig.  2 
also  shows  a  current  path  flowing  upwards  from  the  anode  foil 
through  stack  #2,  and  into  the  cathode  one  jelly  roll  wrap  further 
away.  In  this  work  the  offset  effect  will  be  neglected,  which  is 
reasonable  if  the  number  of  wraps  is  large.  However,  the  double 
stack  effect  must  taken  into  account  as  will  be  described  below. 

Fig.  1  also  attempts  to  show  qualitatively,  using  arrow  length, 
what  would  be  expected  regarding  the  magnitudes  of  the  vari¬ 
ous  currents.  In  particular  the  foil  currents  decrease  in  magnitude 
away  from  the  tabs  as  more  and  more  current  makes  a  right  (or 
left)  turn  and  goes  through  the  stack.  What  is  not  so  obvious 
is  that  the  current  distribution  across  the  stack  is  non-uniform, 
the  current  is  reduced  away  from  the  tabs.  As  will  be  discussed 
below,  this  is  actually  a  transient  effect  at  the  beginning  of  dis¬ 
charge.  The  situation  is  analogous  to  the  current  flow  in  the 
simple  resistor  network  shown  in  Fig.  3.  If  ah  three  vertical  re¬ 
sistors  are  of  equal  value,  and  the  horizontal  resistors  in  any  row 
are  equal,  one  finds  that  less  current  flows  through  the  middle 
resistor  h  <h  =  h- 

3.  Theory 

3.1.  Foil  current  divergence 

Where  possible  units  will  be  written  in  square  brackets,  for 
the  various  symbols  as  they  are  defined.  This  will  hopefully  help 
the  reader  keep  track  of  absolute,  length  and  area  specific  cur¬ 
rents  and  resistances.  Many  expressions  are  duplicated  for  anode 
and  cathode.  In  some  cases  there  are  important  sign  differences 
between  anode  and  cathode,  in  which  case  both  expressions  will 
be  explicitly  written  out  for  clarity. 


Ohms  law  relates  the  current  and  potentials  in  the  foil 

-At/aW  =  -«a/aW  =  v^aW  (1) 

where  for  the  anode  pa  is  foil  resistivity  (Q  cm),  7a  is  the  foil 
thickness  (cm),  Wa  is  the  foil  width  (cm),  0a  is  the  potential 
along  the  foil  (V)  and  R a  =  Pa/TaWa  is  the  length  specific 
foil  resistance  (£2  cm-1),  with  corresponding  definitions  for  the 
cathode.  Note  that  the  possibility  of  distinct  foil  widths  for  anode 
and  cathode  is  included,  this  is  indeed  the  case  in  commercial 
cells. 

One  can  define  the  potential  across  the  stack  at  any  point 
along  the  foils  as 

000  =  0c  00  -  0 aW. 

Similarly  one  can  define  a  stack  current  density  i  (A  cm-2). 
Away  from  the  tabs,  detailed  current  balance  in  the  anode  foil 
gives 

V/aW  =  —Wi(x),  i  /  ij'.  (2) 

where  v-7  is  the  location  of  the  jth  tab,  and  W  =  minimum 
{Wa,  Wq}  is  the  active  width  (usually  the  cathode  is  narrower 
than  the  anode  so  If  =  Wc).  This  just  says  that  away  from  the 
tabs,  the  divergence  of  the  current  in  the  foil  is  proportional  to 
the  flux  across  the  stack  at  that  point.  Right  at  a  tab  the  current 
is  non-analytic,  this  will  be  discussed  in  more  detail  below.  Due 
to  our  choice  of  the  current  direction  in  the  foils,  for  the  cathode 
foil  current  divergence  the  sign  is  reversed 

V/c(x)  =  Wi(x),  i  (3) 


3.2.  Effective  stack  impedance 


In  order  to  proceed  to  the  next  stage,  a  relationship  is  re¬ 
quired  between  the  stack  potential  and  stack  current  density.  A 
comprehensive  description  of  this  relationship  is  provided  by 
the  one-dimensional  stack  models  described  in  ref.  [1].  Because 
of  diffusion  effects  and  active  particle  surface  reaction  kinetics, 
this  relationship  is  in  general  non  linear: 

<p(x)  =  V0(q(x))  -  Rstnck(x)i(x)  +  0(;2)  H - .  (4) 


where  V0(q)  (V)  is  the  open  circuit  voltage  (OCV)  of  the  stack, 
q  (A  h  cm-2)  is  the  state  of  charge  in  the  stack,  Rstack  (^  cm2)  is 
the  effective  stack  impedance  at  i  =  0,  and  0(/2)  indicates  terms 
of  order  i2  or  higher.  V0  can  be  measured  directly  on  a  full  cell 
at  low  current,  or  it  can  be  calculated  from  the  voltage  curves 
of  the  individual  electrodes.  In  the  latter  case  knowledge  of  the 
cell  balance  is  also  required. 

It  is  instructive  to  make  an  Ohmic  approximation  and  assume 
R stack  is  a  constant  in  space  and  time  and  higher  order  terms  in 
(4)  are  neglected.  In  this  case  the  mechanics  of  foil  potentials  and 
currents  are  more  easily  understood.  The  general  case  will  be 
described  in  Appendix  A.  When  higher  order  terms  are  neglected 
one  can  invert  (4) 


V0(q(x ))  -  (p(x ) 


''stack 
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A  rough  estimate  for  the  area  specific  stack  impedance  can  be 
obtained  from  a  volume  weighted  sum  of  electronic  and  ion 
contributions.  For  one  electrode 

^electrode  ~  T 

where  T  is  the  electrode  thickness,  a  is  the  electronic  con¬ 
ductivity  of  the  solid  phase,  k  is  the  conductivity  of  the  elec¬ 
trolyte,  and  the  Us  are  the  corresponding  volume  fractions. 
If  the  volume  fraction  of  inactive  materials  is  negligible  then 
6s°iid  +  ePore  1 .  The  wetted  separator  impedance  is  similar, 
but  without  the  solid  phase  term.  Using  conductivities  from 
[5,6]  for  25  p,m  separator  and  100  p,m  electrodes,  one  finds 
^ stack  ^  2—5  ^2  cm  . 

Eq.  (6)  ignores  the  effective  impedance  of  the  charge  transfer 
reactions  at  the  surface  of  the  active  particles.  Analytic  rela¬ 
tions  for  the  effective  stack  impedance  including  charge  trans¬ 
fer  effects  have  been  reported  elsewhere  [7,3,8].  Using  conduc¬ 
tivity  and  exchange  current  values  from  ref.  [5]  for  a  100  |jim 
thick  electrode,  one  finds  that  the  charge  transfer  resistance 
is  on  the  order  of  3  Q  cm2  per  electrode.  In  total  one  expects 
7? stack  ~  10  £2  cm2  using  this  sort  of  calculation.  Measured  val¬ 
ues  for  power  optimized  cells  [9]  are  in  the  range  30-50  Q  cm2. 

The  7? stack  calculations  described  above  are  for  one  stack 
sandwich.  Since  the  real  cell  has  two  of  these  in  parallel,  it  sees 
half  the  calculated  stack  impedance.  In  addition,  when  inter¬ 
preting  stack  impedances  for  spiral  wound  cells  reported  in  the 
literature,  these  reported  numbers  are  usually  for  two  stacks  in 
parallel. 


where 

9  W(Ra  +  Rc) 

a 2  =  v  (9) 

''stack 

3.4.  Segmented  solution 

Every  unique  tab  position  is  a  point  source  or  sink  for  current. 
The  degenerate  case  where  two  tabs  are  directly  opposite  each 
other  (one  on  the  anode  and  one  on  the  cathode)  amounts  to 
only  one  unique  tab  position.  Tabs  at  the  ends  of  the  electrodes 
(x  =  0  or  v  =  L)  are  referred  to  as  external  tabs.  The  v  interval 
between  any  two  neighboring  tab  positions  will  be  referred  to 
as  a  segment.  For  example,  in  Fig.  1  above,  the  cathode  foil  has 
one  internal  tab  and  no  external  tabs,  the  anode  foil  has  two 
external  tabs  and  no  internal  tabs.  Both  anode  and  cathode  have 
two  segments,  one  from  0  to  L/2,  and  the  other  from  L/2  to  L. 
In  particular  there  is  a  segment  boundary  on  the  anode  directly 
opposite  the  cathode  tab. 

Defining  Lj  =  xj+\  —  xj  as  the  length  of  the  foils  in  the  v 
direction  for  segment  j,  the  dimensionless  group  aLj  describes 
the  in-plane  foil  impedance  relative  to  the  through-plane  stack 
impedance,  and  determines  the  curvature  of  the  foil  current 
profiles  in  segment  j.  When  aLj  <<C  1  the  through-plane  stack 
impedance  dominates  the  cell  impedance,  consequently  the  in¬ 
plane  foil  current  profiles  will  be  linear,  and  the  through-plane 
stack  current  distribution  will  be  close  to  constant.  Conversely 
when  aLj  1  the  in-plane  foil  impedance  dominates,  the  foil 
current  profiles  will  be  curved  and  the  through-plane  stack  cur¬ 
rent  will  be  peaked  at  the  tabs. 


3.3.  Potential  equation 


3.5.  General  solution 


In  general  all  currents  and  potentials  are  functions  of  position 
and  time  (x  and  t).  For  brevity  below  it  is  assumed  this  is  implied 
unless  required  or  otherwise  stated.  Differentiating  (1)  and  using 
(2)  and  doing  similarly  for  the  cathode 

RAWi  =  V20a, 

-RcWi  =  V2</>c. 

Subtracting  gives 

V2</>c  -  V2</>a  =  V2</>  =  -  W(Ra  +  Rc%  (7) 


A  solution  of  (8)  is  called  the  particular  solution,  0p.  Any 
linear  combination  of  the  solutions  to  the  homogeneous  equation 

(V2  —  a2)(p  =  0  (10) 

can  also  be  added  to  0p.  When  a  is  constant,  exp(zbo'v)  are 
two  linearly  independent  solutions  to  (10)  and  the  homogenous 
solutions  will  the  be  same  for  all  segments. 

Anticipating  the  use  of  external  current  boundary  conditions, 
for  calculational  convenience  homogeneous  solutions  are  calcu¬ 
lated  with  unit  slope  at  one  end  of  a  segment  and  zero  slope  at 
the  other  end  of  a  segment,  and  the  particular  solution  with  zero 
slope  at  both  ends.  For  segment  j 


Using  (5)  one  finds 

_ o  Vo  U 

v2d>  =  -W(RA  + 

Tv  stack 

(V2  -a2)<p  =  -a2  Vo 


(V2 

(V2 

(V2 


a2)(p'fl  =  0, 

a2)4>?ht  =  0, 


left 


X=Xj 


=  1, 


left 


x= 


xl+ 1 


V<pfht 


X-Xj 


=  0,  ht 


=  0 


=  1 


a2)^+a2V0  =0,  Vcp] 


=  0,  V(/>] 


X=Xj 


X=Xj+ 1 

=  0. 


*=*;+ 1 

For  constant  a  (Ohmic  condition)  one  finds 
^leit  _  sinh(a[Xy+i  —  x]) 


(8) 


a  sinh(o'[Xj+i  -  xj\) 

right  _  sinh(a[xj  —  x]) 

J  a  sinh(a,|x;-+i  —  xj \ ) 


(11) 
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While  analytic  techniques  exist  for  solving  (11)  for  the  Ohmic 
case  (operator  inversion  methods)  the  resulting  solution,  ex¬ 
pressed  as  a  double  integral,  must  be  evaluated  numerically. 
In  practice  it  is  much  easier  to  solve  (11)  numerically  using  a 
tridiagonal  solver. 

Now  the  general  solution  for  any  segment  can  be  written  as 

<t>i  =  <Pj+  aj(Plft  +  bj(PTjgH 

where  the  constants  aj  and  bj  must  be  determined  by  boundary 
conditions.  Due  to  the  special  properties  of  0ieft  and  0r ight  the 
external  current  coming  in  from  the  left  end  of  the  segment  will 
determine  aj  and  current  from  the  right  will  determine  bj.  Be¬ 
cause  of  the  special  choice  of  boundary  conditions,  the  physical 
interpretation  of  the  particular  solution,  0P,  is  as  follows: 


Since  0c 00 —  0a 00  =  0CO  for  all  v  one  finds  that 
-RaIa(xo)  =  -Rclcixo)  =  cand0A(vo)  =  0cC*o)  =  d.c  has 
dimensions  of  electric  field,  i.e.  (V  cm-1).  Defining  the  relative 
foil  resistances 


Ra 

r  a  =  - 

Ra  +  Rc 

Rc 

rc  =  - 

Ra  +  Rc 
results  in 


(17) 

(18) 


(pA  —  ~  ^a(0P  +  <3</>left  +  bright)  +  CX  +  d  (19) 

(pC  —  rc(0p  +  <J(pMt  +  bright)  +  CX  +  d  (20) 


•  Initially  when  V0  is  constant  then  0P  =  V0. 

•  After  some  current  has  passed  V0  becomes  non  uniform.  Now 
0P  —  Vo  corresponds  to  potentials  (and  therefore  currents) 
along  the  foils  resulting  only  from  the  changes  in  state  of 
charge  along  the  length  of  the  segment,  the  so  called  internal 
mixing  currents  [10].  Mixing  currents  have  the  interesting 
property  that  they  keep  running  after  the  external  current  is 
turned  off. 


3.6.  Electrode  potentials  and  currents 


Since  all  the  physical  boundary  conditions  are  defined  in 
terms  of  the  individual  electrode  potentials,  (0a,  0c)  and  cur¬ 
rents  (I a,  7c)  rather  than  the  stack  potential  and  current,  it  is 
useful  to  derive  expressions  for  the  electrode  potentials  and  cur¬ 
rents.  Using  (7)  and  briefly  dropping  all  segment  indices  one 
finds 

V2<p 


i(x)  = - 

W(RA  +  Rc) 

and  putting  this  into  (2)  and  (3)  and  integrating  one  finds 

1 


(12) 


rx 

IA  =  -W  /  i(x') Ax'  = 

JXo  Ra  +  Rc 

and 

IC  =  W 


V0  +  7a(v0) 


(13) 


'x  1 

i(x)  dx'  = - 

Ra  +  Rc 


V0  +  7c(x0) 


(14) 


where  7A(v0)  and  Iq (jc0)  are  integration  constants.  One  can  now 
plug  these  into  integrated  forms  of  (1) 


<Pa  =  ~Ra  I  Ia(x)Ax 

X0 


Ra 


Ra  +  Rc 


<P  -  RAIA(x0)x  +  (pA(x0 ) 


(15) 


and 

rx 

(pc  =  -Rc  I  lc(x')dxr 

Rc 


X0 


Ra  +  Rc 


0  -  Rcic(xo)y  +  0c(*o) 


(16) 


1  c 

7A  —  ~  ■  ~  (^ 0p  +  a+7 (p\Q ft  +  bV(PYight)  —  (71) 

Ra  +  Rc  Ra 

1  c 

7c  =  ~  '  ~  (^  0p  +  <7^01eft  +  ^V0right)  Xj  (22) 

Ra  +  Rc  Rc 

There  are  now  four  constants  (a,  b,  c,  d)  which  need  to  be  de¬ 
termined  by  external  boundary  conditions.  Each  segment  has  its 
own  profiles  for  potential  and  current,  and  therefore  its  own  set 
of  constants  (a,  b,  c,  d). 

3. 7.  Boundary  conditions 

3.7.1.  Categorizing  tabs 

The  number  of  tabs  for  each  type  are  defined  in  Table  1 .  The 
n1  internal  tabs  are  located  at  positions  j  =  1 ,  ,n1.  The 

electrode  potentials  above  are  only  valid  in  the  regions  between 
these  tabs.  Right  at  the  tabs  0A  and  0c  will  be  continuous  but 
non  analytic,  i.e.  the  first  derivatives,  the  in-plane  foil  currents 
I  A  and  7c,  are  not  always  continuous.  Physically  the  in-plane 
foil  currents  will  exhibit  a  step  change  (usually  changing  sign 
as  indicated  by  the  arrows  in  Fig.  1)  at  each  tab  attached  to 
that  electrode.  The  full  potential  must  be  constructed  from  n1  + 
1  segments,  so  4 (n1  +  1)  boundary  conditions  are  required  to 
determine  all  the  a  j,  bj ,  Cj ,  d/s  (using  j  as  a  segment  index). 
The  problem  is  now  reduced  to  that  of  bookkeeping  for  all  the 
boundary  conditions. 

3.7.2.  Global  boundary  conditions 

For  the  first  segment,  j  =  0,  one  can  arbitrarily  set 

0ao(O)  =  0  (23) 

since  only  potential  differences  are  important.  Alternatively  one 
can  just  set  do  =  0.  Also  the  total  current  flowing  into  (or  out 

Table  1 


Tab  counting 


Tab  type 

Internal 

External 

Total 

Anode 

Cathode 

Both 

nA 

nlc 

n1 

0  <«'  <2 

o  <  4  <  2 

0  <  ne  <  4 

=  +A  +  -  1 
nc  =  nlc  +  nQc  >  1 

n  =  nA  +  nc  =  n1  +  ne 
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of)  each  electrode  is  fixed 

nA 

(lAj(xj)-IAj(xj-i))=-I  (24) 

je anode 
nc 

J2  (Icj(xj)-Icj(xj-i))  =  I.  (25) 

je  cathode 

These  two  are  linearly  dependent,  so  they  only  count  as  one 
boundary  condition.  The  impedance  of  the  tab  material  is  as¬ 
sumed  to  be  negligible  so  the  potentials  of  all  tabs  on  a  given 
electrode  are  all  equal 

4>a(x')  =  <pA(xk),  j  /  k,  j,  k  e  anode 

<Pc(xj)  =  (pc(xk),  j  ±  k,  j,  k  e  cathode. 

There  will  be  n  —  2  such  conditions.  When  there  are  no  external 
tabs  the  currents  at  the  boundaries  will  all  be  0 

/A(0)  =  7C(0)  =  IA(L)  =  IC(L )  =  0. 

When  some  external  tabs  are  present  there  will  be  4  —  ne  such 
conditions. 

3.7.3.  Local  boundary  conditions 

At  each  of  the  tab  positions,  Xj,  the  potentials  (but  not  the 
current)  on  the  tabbed  electrode  must  be  joined.  At  the  corre¬ 
sponding  segment  boundary  on  the  opposing  electrode  (assum¬ 
ing  there  is  no  tab  there)  the  potential  and  current  must  be  joined. 
For  example  if  there  is  an  anode  tab  at  xj  then  both  the  anode 
and  cathode  potentials,  and  the  cathode  current  must  be  joined 
or  matched 

4>Cj-i(xj)  =  (pcj(xj)j  e  cathode 
<(>Aj-\(xj)  =  (pAj(xj)j  e  cathode 
IAj-i(xj)  =  lAj{xj)j  €  cathode 
<t>Aj-\(xj)  =  (pAj(xj)j  e  anode 
4>C j- 1  (x j)  =  4>Cj(xj)j  e  anode 
Icj-i(xj)  =  Icj(xj)j  e  anode. 

There  will  be  3m1  such  conditions.  So  in  total  there  are 
2  +  n— 2  +  4  —  ne  +  3  nl  =  4  +  An1 
boundary  conditions  as  required. 


Cathode 

Tab 

*0  I  *2 


^ -  .segment  u  - ► 

t  J:/  t 


v  Anode  Anode 

Tab  Tab 


Fig.  4.  Schematic  showing  two  segments  and  three  segment  boundaries  3co,  x\ 
and  X2  ■ 

3.7.4.  Example  of  boundary  conditions 

In  order  to  clarify,  all  boundary  conditions  for  the  tab  ar¬ 
rangement  in  Fig.  1  will  be  listed.  To  facilitate  this,  the  two 
segments  and  three  segment  boundaries,  xo,  x\  and  X2,  are 
shown  in  Fig.  4.  All  boundary  conditions  are  listed  in  Table 
2.  Eqs.  (19)— (22)  can  substituted  into  the  relations  in  the  last 
column  of  Table  2,  to  obtain  a  linear  system  of  equations  for 
{<20,  bo,  Co,  do,  a\,  b\,  C\,d\). 

3.7.5.  Equation  system 

Inserting  the  potential  and  current  expressions  (19)-(22)  into 
the  boundary  conditions  results  in  a  system  of  algebraic  equa¬ 
tions  fortheay,  bj,  cj,  dj  coefficients.  This  system  can  be  solved 
manually,  numerically  or  using  symbolic  algebra  system,  de¬ 
pending  on  the  number  of  segments  involved  and  the  desired 
output.  In  this  work  the  numerical  approach  is  preferred  for  the 
following  reasons: 

•  One  must  resort  to  numerical  methods  anyway  for  Li-ion  sys¬ 
tems  where  the  OCV,  V0(q),  is  a  non-trivial  measured  func¬ 
tion,  or  when  the  stack  impedance  becomes  non  trivial. 

•  I  do  not  find  the  messy  analytic  expressions  resulting  from 
solving  a  two  or  more  segment  problem  very  enlightening. 

3.8.  Cell  voltage  and  impedance 

Once  values  for  all  the  aj,  bj,  cj,  dj  coefficients  have  been 
determined,  one  can  graph  and  interpret  internal  quantities,  i.e. 
currents  and  potentials  for  both  the  electrode  foils  and  the  stack. 
These  internal  quantities  are  not  very  amenable  to  experimental 


Table  2 


All  boundary  conditions  for  the  three  tab  cell  design 


Type 

Electrode 

Segment(s) 

Location(s) 

Condition 

End  current 

End  current 

External  current 

Potential  join 

Potential  join 

Current  join 

Tab  potentials 

Arbitrary  potential 

Cathode 

Cathode 

Cathode 

Cathode 

Anode 

Anode 

Anode 

Anode 

0 

1 

0,  1 

0,  1 

0, 1 

0, 1 

0, 1 

0 

xo 

*2 

Xl 

XI 

Xl 

Xl 

XO,  Xl 

X0 

I°c(x  0)  =  0 

Iq(X2)  =  0 
/^Oi)  -  I^(x l)  =  I 

</>c(xi)  =  fc(x\) 
<P°A(x\)  =  f\(x\) 
CA(xi)  =  l\fx  l) 
0a(Xo)  =  0a(*2) 

0a(xo)  =  o 
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verification.  This  means  that  the  calculation  is  the  most  practical 
method  of  accessing  this  information,  but  only  if  it  is  right!  In 
order  to  decide  one  must  calculate  the  external  cell  voltage  and 
impedance  and  compare  with  experiment  for  a  cell  design  for 
which  samples  or  data  are  available. 

The  cell  voltage  is  simply  the  potential  difference  between 
any  anode  tab  and  any  cathode  tab 

Vceii  =  0c (yd  -  0A (yj),  1  e  cathode,  j  e  anode. 


One  can  define  two  distinct  versions  of  cell  impedance,  finite 
difference 


yceii(/  =  0)  -  Ken  CO 

/ 


and  differential 


cdiff  _  dVcell 
Kcell  —  nr 


When  approximating  the  stack  impedance  as  an  Ohmic  or  con¬ 
stant  value,  these  two  versions  of  the  cell  impedance  should 
be  identical.  It  also  turns  out  that  under  these  conditions  the 
cell  impedance  remains  constant  during  the  discharge.  The  fi¬ 
nite  difference  version  is  fairly  straightforward  to  calculate  but 
requires  one  to  keep  track  of  the  7  =  0  solutions  to  all  the  equa¬ 
tions  above.  Calculating  the  differential  form  requires  the  so¬ 
lution  of  a  second  system  of  boundary  condition  equations  that 
are  implicitly  differentiated  with  respect  to  7.  The  solution  yields 
ddj/dl,  dbj/dl,  dcj/dl,  ddj/dl. 


3.9.  Foil  heating 


In  high  power  applications  it  is  essential  to  understand  heat 
generation  effects  at  high  currents.  Total  heat  is  easily  calculated 
from  I2  R.  However,  lot  of  information  about  the  cell  internals 
is  also  available,  which  can  be  used  to  learn  where  the  heat  is 
generated.  This  is  another  effect  that  is  very  difficult  to  measure 
directly.  The  local  heating  rate  in  the  anode  foil  will  be 
1  ? 

?aW  =  —IaR  a 
WA 

with  a  similar  expression  for  the  cathode.  qA  will  have  dimen¬ 
sions  of  (W  cm-2).  The  local  heating  rate  in  the  stack  is 

^/stackC^)  =  l  ^ stack 

with  the  same  dimensions  as  The  local  heating  rates  can 
easily  be  integrated  over  space  to  get  total  heating  rates,  and/or 
over  time  to  get  total  heat  generated.  In  order  to  keep  things  as 
simple  as  possible,  reversible  or  entropic  heating  (or  cooling)  is 
not  included  in  the  stack  heat  generation. 

It  is  worth  pointing  out  that  since  the  foil  currents  are  usually 
peaked  near  the  tab  positions,  the  heat  generation  rate  in  the  foils 
will  be  even  more  strongly  peaked  near  the  tabs.  For  example 
even  when  aL  <<C  1  and  the  foil  currents  are  linear,  the  heat 
generation  will  decay  quadratically  moving  away  from  a  tab. 
This  will  be  shown  graphically  below. 

3.10.  Time  dependence 


All  of  the  above  theory  describes  the  current  and  potential 
distributions  at  the  instant  the  current  is  turned  on.  As  the  cell 


discharges  all  internal  quantities  acquire  a  time  dependence,  and 
the  state  of  charge  in  the  stack  varies  according  to 

q(x,t)=  f  i(x,t')dt' 

Jo 

which  quickly  induces  non-uniformity  in  V0(q).  Non- trivial  ef¬ 
fects  can  arise  from  the  shape  of  V0(q).  When  it  is  flat  as  a 
function  of  q  (voltage  plateau)  then  large  composition  gradients 
can  develop  along  the  length  of  the  cell  (x  direction).  When  V0(q) 
is  steeply  sloping,  any  existing  composition  gradients  along  the 
length  will  be  quickly  eliminated  by  internal  mixing  currents 
arising  from  0p. 

In  practice  a  numerical  simulation  is  required  using  stan¬ 
dard  finite  element  methods.  Results  of  such  simulations  for 
an  LixMn204/Lii_xC6,  spiral  wound  Li-ion  cell  will  be  de¬ 
scribed  in  the  next  section.  The  goal  is  to  demonstrate  some 
of  the  interesting  effects  that  can  occur  with  the  foil  currents. 
For  this  purpose  it  is  best  to  use  a  constant  value  for  the  stack 
impedance. 

4.  Example  calculations 

4.1.  Cell  details 

Results  will  the  shown  for  two  cell  designs  using  the  non 
trivial  tab  arrangement  shown  in  Fig.  1 .  Cell  #1  is  intended  to  be 
towards  the  power  optimized  end  of  the  spectrum,  and  cell  #2 
is  intended  to  be  energy  optimized.  Design  parameters  for  both 
cells  are  shown  in  Table  3. 

The  foil  resistivities  p\  and  pc  correspond  to  Cu  and  Al  re¬ 
spectively.  The  electrode  widths  are  typical  of  those  found  in 
commercial  18650  size  Li-ion  cells.  The  power  optimized  cell 
(cell  #1)  uses  thicker/longer  foils  and  a  thinner  stack,  which 
is  indirectly  expressed  through  the  smaller  stack  resistivity.  As 
discussed  above  the  larger  value  for  aL  in  cell  #1  results  in  a 
non-linear  shape  for  the  foil  currents.  In  order  to  get  some  sig¬ 
nificant  polarization  along  the  electrode  foils,  a  sizeable  current 
must  be  applied.  In  the  examples  below  the  current  is  50  A  of 
discharge  current. 

The  OCV  curve  was  calculated  by  differencing  analytic  ap¬ 
proximations  to  V(q)  (described  in  the  appendix  in  ref.  [11]) 
for  the  anode  and  cathode  active  materials.  The  grid  spacings 
Ax  =  0.1  cm,  At  =  0.5  s  were  used.  Because  there  are  no  di¬ 
rect  comparisons  to  a  real  cell  in  this  work,  the  results  for  V^n 

Table  3 


Example  cell  designs 


Parameter 

Units 

Cell  #1 

Cell  #2 

cm 

5.7 

5.7 

7a 

pm 

15 

10 

Pa 

p£2  cm 

1.7 

1.7 

Ra 

m£2  cm-1 

0.20 

0.30 

WC 

cm 

5.5 

5.5 

7c 

pm 

20 

15 

Pc 

p£2  cm 

3.1 

3.1 

Rc 

m£2  cm-1 

0.28 

0.38 

L 

cm 

200 

100 

A  tack 

£2  cm2 

5 

50 

olL 

— 

4.60 

0.86 
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Fig.  5.  Cell  #1  foil  potentials  0  a (jc)  and  0c (x). 


and  RCq\\  will  not  be  discussed.  Also  when  Rsta ck  is  constant 
their  behavior  is  trivial. 

4.2.  Instantaneous  results 

In  this  section  results  are  shown  at  t  =  0,  the  instant  the  cur¬ 
rent  is  turned  on.  At  this  point  in  time  the  electrons  are  assumed 
to  rearrange  to  their  equilibrium  charge  density  distribution,  but 
no  ions  have  moved.  The  calculated  foil  potentials  for  cell  #1 
are  shown  in  Fig.  5.  In  all  such  figures,  the  vertical  dashed  line  at 
x  =  L/  2  indicates  the  position  of  the  cathode  tab.  It  should  be  re¬ 
membered  that  anode  tabs  are  located  at  the  left  and  right  edges, 
v  =  0  and  x  =  L.  The  potentials  for  cell  #2  are  qualitatively 
similar  and  are  therefore  not  shown.  The  anode  foil  potential  is 
highest  at  the  tabs  and  exhibits  a  minimum  half  way  between 
the  two  tabs.  The  cathode  foil  potential  is  a  minimum  at  the  tab. 

The  foil  currents  for  both  cells  are  shown  in  Figs  6  and  7. 
When  viewing  these  graphs  keep  in  mind  that  positive  currents 
are  moving  to  the  right  in  Fig  1  and  negative  are  moving  to  the 
left.  Because  of  the  symmetry  in  the  cell  designs  considered,  the 
current  entering  each  anode  tab  is  exactly  1/2  =  25  A.  In  the  case 
of  broken  symmetry,  i.e.  moving  the  cathode  tab  off  center,  the 
theory  would  tell  us  exactly  how  the  current  would  be  unevenly 
distributed  between  the  two  anode  tabs.  At  the  cathode  tab  one 
sees  a  discontinuity  in  the  current,  jumping  from  25  A  to  —25  A. 
For  internal  tabs  such  current  discontinuities  are  more  the  rule 
than  the  exception.  One  can  readily  see  that  the  foil  currents  for 
cells  1  and  2  are  qualitatively  different.  For  cell  #2  the  currents 
show  linear  decay,  as  expected  since  aL  is  small,  whereas  the 
currents  for  cell  #1  with  larger  aL,  are  noticeably  non-linear.  For 
cell  #1  one  can  also  see  some  asymmetry  within  each  segment, 
the  currents  decay  faster  near  the  cathode  tab  because  Rq  >  Ra- 
This  asymmetry  is  also  evident  in  the  stack  current  shown  in  Fig. 
8  which  is  more  peaked  near  the  cathode  tab  than  at  either  anode 
tab. 

4.3.  Dynamic  effects 

As  alluded  to  above  the  high  stack  current  at  the  tab  locations 
seen  in  Fig.  8,  is  a  transitory  situation.  As  the  discharge  proceeds 


Fig.  6.  Cell  #1  foil  currents  7a  (x)  and  7c  (x). 


Fig.  7.  Cell  #2  foil  currents  7a  (x)  and  7c  (x). 


the  state  of  charge  (SOC)  becomes  non-uniform  as  a  function 
of  v.  This  in  turn  results  in  a  non-uniformity  in  the  open  circuit 
potential  V0(q(x)),  with  higher  voltages  away  from  the  tabs.  The 
voltage  effect  will  induce  in  increase  in  the  stack  current  away 
from  the  tabs.  Finally,  near  the  end  of  discharge  the  active  ma¬ 
terial  near  the  tabs  will  become  used  up  (discharged),  at  which 
time  all  current  must  flow  through  the  stack  away  from  the  tabs 
where  some  capacity  remains.  All  of  these  effects  are  shown  in 


Fig.  8.  Cell  #1  stack  current  i(x). 


J.N.  Reimers  /  Journal  of  Power  Sources  158  (2006)  663-672 


671 


Fig.  9.  Cell  #1  dynamic  evolution  of  the  stack  current,  i(x,  t). 


Length  (cm) 

Fig.  10.  Cell  #1  dynamic  evolution  of  the  foil  currents,  /a(v,  t)  and  /c(v,  t). 

Fig.  9.  Some  of  the  fine  structure  in  these  profiles  is  due  to  the 
non  trivial  shape  of  V0(q). 

Similar  effects  are  also  observed  in  the  foil  currents  (Fig.  10). 
In  particular,  at  10%  SOC  one  can  see  the  cathode  foil  running 
almost  full  current  (25  A)  for  about  40  cm  on  either  side  of  the 
cathode  tab,  before  linear  decay  begins. 


Fig.  11.  Cell  #1  heating  generation  breakdown. 


Fig.  12.  Cell  #2  heating  generation  breakdown. 


Fig.  13.  Cell  #1  dynamic  evolution  of  the  total  heat  generation. 


4.4.  Heating  effects 

A  break  down  of  various  heating  effects  for  both  cells,  at 
t  =  0,  are  shown  in  Figs.  11  and  12.  As  anticipated  significant 
heat  generation  is  predicted  near  the  tabs  in  cell  #1,  to  which 
both  the  stack  and  foil  contribute.  Cell  #1  is  dominated  by  foil 
heat  generation  and  cell  #2  is  dominated  by  more  uniform  stack 
heat  generation  (note  the  different  y-axis  scales).  Fig.  13  gives 
a  flavor  for  the  time  dependence  of  the  total  heat  generation 
for  the  cell  #1.  One  can  see  that  near  the  end  of  discharge  the 
heat  generation  tends  to  move  away  from  the  tabs,  but  one  never 
observes  a  peak  in  the  heating  rate  away  from  the  tabs.  This  is 
just  a  consequence  of  the  dominating  effect  of  foil  heating  in 
this  cell. 

5.  Conclusions 

The  purpose  of  this  work  is  two-fold 

(1)  Outline  the  theory  for  calculating  internal  foil  currents  and 
potentials  for  a  general  cell  design  with  an  unlimited  number 
of  tabs,  randomly  located. 

(2)  Provide  some  qualitative  understanding  of  how  potentials, 
current  and  heating  behaves  during  discharge  of  a  cell. 
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In  cells  with  low  stack  impedance  and  long  electrodes,  polar¬ 
ization  effects  in  the  foils  cannot  be  neglected.  The  value  of  the 
dimensionless  group  aL  is  easily  calculated,  and  can  be  used  to 
decide  if  foil  polarization  is  important  for  any  given  cell  design. 
The  degree  of  cell  heating  caused  by  the  foils  will  in  general 
depend  on  aL,  if  aL  1  the  foil  heating  will  make  a  major 
contribution,  especially  near  the  tabs.  This  heating  will  in  turn 
change  the  kinetics  in  the  stack  through  its  dependence  on  elec¬ 
trolyte  conductivity,  salt  diffusion,  surface  reaction  kinetics  and 
solid  state  diffusion.  Taking  proper  account  of  all  these  inter¬ 
acting  effects  would  require  the  incorporation  of  a  detailed  dy¬ 
namic  model  for  the  stack  impedance.  This  may  be  the  subject 
of  a  future  paper.  I  hope  that  these  ideas  and  methods  will  be  of 
assistance  in  designing  high  power  cells. 

Appendix  A.  Non  Ohmic  stack  impedance 


When  the  relation  between  the  stack  potential  and  stack  cur¬ 
rent  density  is  non-linear,  then  iterative  methods  will  be  required. 
This  is  always  the  case  when  a  detailed  stack  model  is  used  as  de¬ 
scribed  in  ref.  [1].  Because  the  state  of  charge  quickly  becomes 
non-uniform  along  the  length  of  the  electrodes,  one  must  im¬ 
plement  numerous  instances  of  the  full  stack  model  each  main¬ 
taining  its  own  state  information.  The  result  of  this  is  that  the 
stack  impedance  now  becomes  non-uniform  and  time  depen¬ 
dent:  ^  stack  (A?  t).  All  the  boundary  conditions  equations  in  (1 1) 
must  be  solved  numerically  using  a  tridiagonal  solver. 

Because  i  is  not  a  small  parameter,  the  expansion  in  (4)  is  not 
convergent.  Given  the  stack  current  from  the  previous  time  step, 
or  some  initial  guess  iQ  one  can  expand 


(0  —  0o )  +  O(S02)  ~  i0 


1 


^  stack 


(0  -  0o) 


(26) 


This  corrected  stack  current  expression  can  then  be  used  in  (7). 
The  foil  and  stack  model  must  be  iterated  until  self  consistency 
is  achieved.  This  amounts  to  a  decoupled  two-dimensional  sim¬ 
ulation.  Results  of  such  simulations  may  be  presented  in  a  future 
publication. 
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